Weak decay constant of pseudoscalar mesons in a QCD-inspired model* 
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We show that a linear scaling between the weak decay constants of pseudoscalars and the vector 
meson masses is supported by the available experimental data. The decay constants scale as fmj fit = 
Mv/Mp (f m decay constant and Mv vector meson ground state mass). This simple form is justified 
within a renormalized light-front QCD-inspired model for quark-antiquark bound states. 



I. INTRODUCTION 



Effective theories used to describe hadrons, which are 
inspired by Quantum Chromodynamicspl 0, Q can be 
useful in indicating direct correlations between observ- 
ablcs of different hadrons. Therefore, it is possible to pin 
down the relevant dependence of the observables with 
some physical scales that otherwise would have no simple 
reason to present a direct relation, besides being proper- 
ties of the same underlying theory. For example, a sys- 
tematic dependence of a hadron observable with its mass 
can offer useful guide for presenting results obtained in 
Lattice QCD. Systematic correlations between different 
meson properties with mass scales are also found from 
the solution of Dyson-Schwinger equations 4] . 

One intriguing aspect is the dependence of the weak 
decay constant of the pseudoscalar meson with its mass. 
For light mesons up to D, the weak decay constant tends 
to increase with the mass, while numerical simulations of 
quenched lattice-QCD indicate that fry > £biM> which is 
still maintained with two flavor sea quarks yjug. General 
arguments, within Dyson-Schwinger formalism for QCD 
in the heavy quark limit, says that the weak decay con- 
stant should be inversely proportional to V^m0 (Af m 
is the pseudoscalar mass). Effective QCD inspired mod- 
els valid for low energy scales can also be called to help 
to investigate this subtle point. In these models 0,0,01, 
the interaction is flavor independent, while the masses 
of constituent quarks can be changed, which naturally 
implies in correlations between observables and masses. 

Our aim here, is to investigate the pseudoscalar weak 
decay constant within a QCD inspired model @- The ef- 
fective mass operator equation for the lowest Light-Front 
Fock-state component of a bound system of a constituent 
quark and antiquark of masses mi and ni2, obtained in 
the effective one-gluon-exchange interaction approxima- 



tion and simplified in the fj.-model0 0] to 



M^Vm(£,fc_L) = 
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where the phase space factor is 



A ) lfjm(x', k'j_), (1) 
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and ip m is the projection of the light-front wave-function 
in the quark-antiquark Fock-state. The mean square mo- 
mentum transfer [(k[ — k\) 2 + {k' 2 — fc 2 ) 2 )/2 gives — Q 2 
(ki and k[ are the quark four- momenta) . The coupling 
constant a defines the strength of the Coulomb-like po- 
tential and A is the bare coupling constant of the Dirac- 
delta hyperfine interaction. The energy transfer in Q 2 is 
left out. Confinement comes through the binding of the 
constituents in the meson, which in practice keeps the 
quarks inside the mesons. 

The mass operator equation JQ| needs to be regularized 
and renormalized in order to give physical results, such 
development has been performed in Ref . Q . In that work, 
it was obtained the renormalized form of the equation for 
the bound state mass, which is i) invariant under renor- 
malization group transformations, ii) the physical input 
is given by the pion mass and radius Hi) no regularization 
parameter. 

In the work of Ref. |a| , the quark mass was changed to 
allow the study of mesons with one light antiquark plus 
a strange, charm or bottom quark. The masses of the 
constituent quarks were within the range of 300 up to 
5000 MeV. For the up-down quarks a mass of 384 MeV 
is found from the rho meson mass, which in the model is 
weakly bound. The Dirac-delta interaction comes from 
an effective hyperfine interaction which splits the pseudo- 
scalar and vector meson states. In the singlet channel the 
hyperfine interaction is attractive, which is not valid for 
the spin one mesons. In the model, the Dirac-delta in- 
teraction mock up short-range physics which are brought 
by the empirical value of the pion mass, and a reason- 
able description of the binding energies of the constituent 
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quarks forming the pseudoscalar mesons is found 8]. The 
model, without the Coulomb like interaction, was also 
able to describe the binding energies of the ground state 
of spin 1/2 baryons containing two light quarks and a 
heavy oneQ- 

Within the effective model of Eq.Q, the low- lying 
vector mesons are weakly bound systems of constituent 
uarks while the pseudo-scalars are more strongly bound 
This allows to calculate the masses of the constituent 
quarks directly from the masses of the vector mesons 
ground states @: 



1 



m u 


= -M B 
2 p 


= 384 MeV , 




m s 


= M K , 


- -M p = 508 MeV , 




m c 


= M D , 


- -M p = 1623 MeV , 




m b 


= M B , 


- -M D = 4941 MeV , 
2 p 


(2) 



where it is used the values of 768 MeV, 892 MeV, 2007 
MeV and 5325 MeV for the p, K*, D* and B* masses, 
respectively [Tflj. 

Here, we use the effective model to predict a physi- 
cal property directly related to the wave-function of the 
ground state of the pseudo scalar mesons. We calcu- 
late the weak decay constants (/ m ) of K + , D + , Df, for 
which experimental values are known [lOj. Besides the 
constituent quark masses from Eq.Q and the pion mass, 
our calculation needs as input the pion weak decay con- 
stant, f n = 92.4±0.07±0.25MeV 10]. The eigenfunction 
of the interacting mass squared operator from Eq. for 
large transverse momentum behaves as the asymptotic 
wave-function, which decreases slowly as pj 2 . There- 
fore, in the calculation of the weak decay constants it is 
necessary to regulate the logarithmic divergence in the 
transverse momentum integration and take care of the 
cut-off dependence to be able to give an unique answer. 
One has to consider that the pion decay constant pro- 
vides the short-range information contained in the pion 
wave function, which we suppose to be the same for all 
pseudo-scalars. Here, we just write the divergent integral 
in the transverse momentum in terms of f v and from that 
obtain the other decay constants. 



II. MESON LIGHT-FRONT WAVE FUNCTION 

The wave function of the meson (m) is the solution 
of Eq.Q. In the approximation where the Coulomb- 
like interaction is considered in lowest order, the pseudo- 
scalar meson wave function is given byQ 



ip m (x,k±) 



1 - 



1 



Gr, 



Ml 



y/ X (l - x) Ml 

dx'dk'J(x')6{\ - x') / 4mim 2 a 
^/x'(l~x') V 3^ 2 



Ml - M'i 



where 



P 



k\ 



1 



(3) 



(4) 



in the frame in which the meson has zero transverse mo- 
mentum. (M'q is obtained from Mq by substitution of 

k± and x by k' ± and x' , respectively.) The overall nor- 
malization of the qq Fock-component of the meson wave- 
function J2J is G m . 

In this first calculation of the decay constant within 
this model, we are going to assume the dominance of the 
asymptotic form of the meson wave function and simply 
use 

^m(^fcx) = . = = „ 2 G " 2 ■ (5) 



y/x(l - x) Ml 



Mi 



To obtain the pseudoscalar decay constants, we follow 
Ref. [ll| . To construct the observables in terms of the me- 
son wave function, one has to account for the coupling 
of the quark spins, which is described by an effective La- 
grangian density with a pseudo-scalar coupling between 
the quark (qi(x) and q2{x)) and meson ($ m (a;)) fields [Tl| 

-iG m $ m (x) q l (x)-i b q 2 {x 



£eff(x) 



h.c. 



(G) 



the coupling constant is G m . From the effective La- 
grangian above one can derive meson observables and 
write them in terms of the light-front asymptotic wave 
function, Eq. 0. To achieve this goal, it is necessary 
to eliminate the relative x + -timc [x + = t + z) between 
the constituents in the physical amplitude, which then al- 
lows to write the meson observable in terms of the wave 
function 



III. RESULTS FOR THE WEAK DECAY 
CONSTANT OF PSEUDOSCALAR MESONS 

The pseudoscalar meson weak decay constant is calcu- 
lated from the matrix element of the axial current A 11 (0) , 
between the vacuum state |0) and the meson state \q„ 
with four momentum q m [lfjj : 



(0 | A"(0) | q r , 



iV2f„ 



(7) 



where A^(x) —q^-f^^q^). 

Using the pseudoscalar Lagrangian, Eq. 
calculate the matrix element of the axial current, which 
is expressed by a one-loop diagram and can be written 
as: 



iV2M m f„ 



N r G r . 



d A k 

where 7 + = 7 + 7 
Siijp) = i/(jf> - mi + 
field. 



-Tr [7 + 7 5 S 2 (fe)7 5 Si (k- q m )] , (8) 

, N c = 3 is number of colors and 
ie) is the propagator of the quark 
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TABLE I: Results for the pseudoscalar meson weak decay 
constants f m calculated with Eg. 11211 . The inputs for the 
model are m q , m-q, M m and / w given in the table. All masses 
and decay constants are in MeV. ( ' The experimental mass 
of Dg + is quite near to the model cs vector meson mass which 
is given by m c + m s = 2131 MeV. ) 



qq 




M V \W] 


rmodel 
Jm 


/-mo] 


n + (ud) 


140 


771 (p) 


92.4 


92.4 ± .07 ±0.25 


K + (us) 


494 


892 (AT*) 


107 


113.0 ± 1.0 ±0.31 


D+(cd) 


1869 


2010 (D*+) 


241 


91 9 +127+6 
^ lz -106-28 


D+(cs) 


1969 


2112+ (D*+) 


253 


201 ± 13 ±28 



By integration over k~ in Eq. JSJ), the relative light- 
front time between the quarks is eliminated and one ob- 
tain the expression of f m suitable for the introduction of 
the meson light front wave-function. So, performing the 
Dirac algebra and integrating analytically over k~ , one 
obtains 



fn 



V2 
8tt 3 



Nr. 



dki 



dx ((1 — x)m,2 + xmi) 



x(l — x)M^ — k\ — m\x — m2(l — x) 



(9) 



in the meson rest-frame. We have used the momentum 
fraction x = k + 

One can write Eq.@ in terms of the valence compo- 
nent of the pseudoscalar meson wave function as: 



V2 



dx 



fm — „ n N c _ 

8tH J y/ x (l-x) 
x / dk\tp m (x, fcjj ■ 



((1 — x)rri2 ± xm{) 



(10) 



The above expression is general and one can use it to 
calculate the decay constant of any pseudoscalar meson 
state, and as well one can use it to normalize the eigen- 
function of the squared mass operator from the solution 
ofEq.lU). 

We observe that Eq.©, written in terms of the asymp- 
totic part of the valence wave function has a logarithmic 



divergence in the transverse momentum integration due 
to the slow decrease of the wave function. From a physi- 
cal point of view, one could think that the regularization 
scale is larger than the masses of the quarks and the 
divergent transverse momentum integration, will be de- 
fined through the value of f^, for example. Therefore, 
one has: 



f m = const. / dx ((1 — x)m2 + xm\) , (11) 
Jo 



and const, determined by /„•. One observe as well that, 
f m cx mi + m 2 , which in our model is the vector meson 
mass, thus one immediately gets: 



fm 



Mv 

Mn 



(12) 



The numerical results of Ea. (|12J) are shown in Table 
I. It is verified that a reasonable description of the weak 
decay constants of the pseudoscalar mesons is possible 
within the effective light-front model. However, we have 
made use only of the asymptotic form of the wave func- 
tion and one needs to investigate the decay constant with 
more refined wave functions, eigenstates of the squared 
mass operator, Eq.Q, which includes the dynamics of 
the effective quarks. Therefore, the results which are 
overestimating the heavier meson decay constants, can 
be an indication that a more elaborated wave function is 
needed, although one cannot discard that other mecha- 
nisms could be relevant 0- 

Also, we intend to perform the evaluation of the weak 
decay constants using a more sophisticated version of the 
model, where confinement is included |l3j| . which so far 
was shown to describe the meson spectrum. 

In summary, we have shown the existence of a direct 
proportionality between the weak decay constants and 
the masses of the vector mesons ground states, which 
can be an useful tool in the systematic study of these 
quantities. 
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